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Abstract 

The type III Hermite X m exceptional orthogonal polynomial family is generalized 
to a double-indexed one X mi . m2 (with m\ even and rri2 odd such that m% > mi) and 
the corresponding rational extensions of the harmonic oscillator are constructed by 
using second-order super symmetric quantum mechanics. The new polynomials are 
proved to be expressible in terms of mixed products of Hermite and pseudo-Hermite 
ones, while some of the associated potentials are linked with rational solutions of the 
Painleve IV equation. A novel set of ladder operators for the extended oscillators is 
also built and shown to satisfy a polynomial Heisenberg algebra of order mi — m\ + 1, 
which may alternatively be interpreted in terms of a special type of {mi — m\ + 2)th- 
order shape invariance property. 
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1 Introduction 



During the last few years, a lot of research activity has been devoted to the construction 
of exceptional orthogonal polynomials (EOP), which are new complete and orthogonal 
polynomial systems extending the classical families of Hermite, Laguerre and Jacobi. In 
contrast with the latter, the former admit some gaps in the sequence of their degrees, the 
total number of them being referred to as the codimension. 

The first examples of EOP, the so-called Laguerre and Jacobi X\ families, with codi- 
mension m — 1, were proposed in the context of Sturm-Liouville theory [TJ [2]. Shortly 
thereafter, such EOP were applied to quantum mechanics and proved to be related to some 
exactly solvable (ES) rational extensions of well-known potentials [3] . These new quantum 
potentials had the additional interesting property of enlarging the known list of transla- 
tionally shape-invariant (SI) potentials, which was thought to be complete at that time 

In agreement with some previous works on algebraic deformations of SI potentials [SIE], 
it then appeared convenient to use a supersymmetric quantum mechanical (SUSYQM) 
method to construct some additional examples of quantum potentials related to X\ EOP 
[8] and to extend such an analysis by proposing some new Laguerre and Jacobi X2 families 
[2]. The transition from codimension m = 1 to codimension m = 2 was characterized by a 
growing complexity, since there appeared three distinct families of EOP (labelled as types 
I, II and III, respectively), the first two being associated with SI potentials and strictly 
isospectral to the partner conventional potential. 

The next step was the obtention of type I and type II Laguerre and Jacobi X m EOP 
families and related potentials for arbitrarily large codimension m [T0| [TTj and the study 
of their properties [T21 It then became clear that, in the m = 1 special case, the 
polynomials of the type I and type II families happened to be proportional, hence leaving 
a single independent one. In addition, the X m EOP where shown to be obtainable through 
several equivalent approaches, such as the Darboux-Crum transformation [HI [151 HB], the 
Darboux-Backlund one [T7] and the prepotential method |18j . 
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Type III X m EOP families were also studied in detail and shown to exist only for even 
codimension [17]. Some results were also obtained for type III Jacobi ones [IB]. 

The development of these type III EOP families helped in understanding the absence of 
Hermite X% EOP families in the first works carried out on the subject [U(2]- Some Hermite- 
type polynomials had indeed been obtained in quantum mechanical models in the early 90s 
[T91 I2"U] and re-discovered several times later on (see, e.g., [61 EB E2J E3]), where they 
appeared in the vast domain of anharmonic oscillators constructed by SUSYQM techniques 
(see, e.g., [2U [25j [261 [28] and references quoted therein). It became clear [291 [30] that 
such Hermite-type polynomials actually belonged to type III EOP families corresponding 
to even codimension and that no other Hermite EOP families could be constructed in a 
first-order SUSYQM approach (nor in any of the above-mentioned equivalent methods). 

Recently, it was shown that the list of X m EOP families (and associated potentials) was 
not exhaustive and that novel multi-indexed ^ mi ,m 2 ,...,m fc families could be obtained through 
the use of multi-step Darboux algebraic transformations [31], the Crum-Adler mechanism 
[32J , higher-order SUSYQM [331 EI] or multi-step Darboux-Backlund transformations [35] . 
These works aimed at generalizing type I and type II X m EOP families and consequently 
at extending the class of SI potentials that are rational extensions of conventional ones. 

It is worth mentioning too that very recently some studies were devoted to the construc- 
tion of new rational extensions of solvable potentials with a finite bound-state spectrum 
and to the study of the polynomials appearing in their wavefunctions [291 EQl ESI E3 138] . 
as well as to some mathematical aspects of the theory [39l 140] . 

As far as the authors know, however, no attempt has been made so far at extending 
type III X m EOP families to multi-indexed ones. The purpose of the present paper is to 
start filling in the gap by considering the simplest example, namely Hermite X mii?n2 EOP 
families that may occur in a second-order SUSYQM approach to the harmonic oscillator. 
Although the latter subject has a long history (see, e.g., [261 123 [28] and references quoted 
therein), the EOP viewpoint will shed some new and interesting light on it. 

Another purpose of this paper is to re-examine the construction of ladder operators 
for harmonic oscillator partners. We plan to show that for the potentials associated with 
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Hermite X mij?n2 EOP, there exists an alternative of the usual procedure that combines 
the oscillator creation and annihilation operators with the supercharges [2U ESI EZ] or 
combinations of supercharges [28j ED H2]- 

In section 2, we review the construction of Hermite X m EOP families and associated 
harmonic oscillator rational extensions in first-order SUSYQM. We extend both of them to 
second-order SUSYQM in section 3. Ladder operators are then built for the corresponding 
potentials in section 4. Finally, section 5 contains the conclusion. 

2 Harmonic oscillator rational extensions in first- 
order SUSYQM 

In first-order SUSYQM [5J, one considers a pair of SUSYQM partners H^, 

if« = A* A = — f-, + V {+ \x) - E, ifH = AA^ = — f-, + V { '\x) - E, 
dx 2 ax 2 

A^ = — < L + W(x), A = 4- + W(x), t 2 " 1 ) 
ax ax 

V {±) (x) = W 2 (x) =F W'(x) + E, 

which intertwine with the first-order differential operators A and A^ as AH^ = H^A and 
A^H^ = H^A^. Here W(x) is the superpotential, which can be expressed as W(x) = 
— (log <j)(x))' in terms of a nodeless seed solution <p(x) of the initial Schrodinger equation 

+ U (+) (x)U(x) = Ecf)(x), (2.2) 
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dx 

E is the factorization energy, assumed smaller than or equal to the ground-state energy of 
V {+ \x), and a prime denotes a derivative with respect to x. 

Whenever both energies are equal and 4>(x) is therefore the ground-state wavefunction 
of V( + \x), the partner potential V^(x) has the same bound-state spectrum as V^ + \x), 
except for the ground-state energy which is removed (case i). For E less than the ground- 
state energy, in which case <p(x) is a nonnormalizable eigenfunction of V^ + \x), the partner 
V^~\x) has the same spectrum as V^ + \x) if <p~ l (x) is also nonnormalizable (case ii) or 
has an extra bound-state energy E below the spectrum of V^ + '(x) corresponding to the 
wavefunction _1 (a;), if the latter is normalizable (case Hi). 



As well known, for the harmonic oscillator potential V(x) = x 2 (x EM), the correspond- 
ing Schrddinger equation has an infinite number of bound-state wavefunctions, which can 
be written as 



i/; u (x) =KH u (x)e ^ 



M v = (^2 v u\)~ 1/2 , u = 0,1,2, 



(2.3) 



where H v (x) is a z/th-degree Hermite polynomial [43J. The associated bound-state energies 

are given by 2v + 1, v = 0, 1, 2, 

If one assumes 

V( + \x)=x 2 (2.4) 

and E = 1 in (12. ip and (12. 2p . then W(x) = x and the partner potential V^~'(x) = x 2 + 2 
is a translated oscillator, which reflects the SI of the harmonic oscillator [5] . 

If now E < 1, the only possible polynomial- type nodeless seed solutions of (12 .2p are 
provided by the functions 



(x) = H m (x)e^ x2 , 



m = 2,4,6, 



(2.5) 



where T-L m (x) is a mth-degree pseudo-Hermite polynomial [43 

[m/2] 



defined by 



H m (x) = {-i) m H m (ix) = ml u 1 o M 

L — ' p\[m — 2p)\ 



m—2p 



p=0 



(2.6) 



The corresponding factorization energies are 



E„ 



-2m - 1. 



(2.7) 



Since 0~ x (x) is normalizable on M, we are in case in of SUSYQM. It is worth observing 
here that for odd m, the function (f) m (x), as defined in (12.51) and (I2.6p . is also a solution of 
equation (12.21) corresponding to (12.71) . but as it vanishes at x — 0, it does not qualify as a 
seed function in first-order SUSYQM. 

The corresponding partner potential is now 



V { - ] (x) = x 2 - 2 



nil, 
T-ir 



nil 

rL rr, 



(2.8) 



^^This pseudo-Hermite polynomial may also be seen as a special case of generalized Hermite polynomial 

SI ESI. 
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and the spectra of and ) are 

El +) =2(z/ + m + l), 1/ = 0,1,2,..., (2.9) 

and 

BH=2(i/ + m + l), i/ = -m- 1,0,1,2,..., (2.10) 

respectively. 

The excited-state wavefunctions of if (~- ) can be found by acting with the operator A, 
where 

W(x) = -x-^, (2.11) 

on the wavefunctions ipi + \x) = if) u (x) of ffW. On using the Hermite and pseudo-Hermite 
polynomial identities given in appendix A, we obtain 

4~\x) = MH^^J™l +1 ( x ), v = 0, 1, 2, . . . , 

yi+L+i^) = -^(xjff.+iix) - 2mU m -i(x)H l/ (x), (2.12) 

atw = ^= = ^ f = = ytr+\ v + m + m- l/2 . 

J E (+) V 2 ( z/ + m + 1 ) 



On the other hand, the ground-state wavefunction can be written as [23] 

-^J . (2.13) 

Hence the whole set of wavefunctions is given by 

^-\x) = J^-)±-L-y^\ x ), n = u + m + l, v = -m - 1, 0, 1, 2, . . . , (2.14) 
U m [x) 

where we define 

yt\x) = 1. (2.15) 

The set of nth-degree polynomials yn(x), n — 0, m + 1, m + 2, . . . , being orthogonal 
and complete with respect to the positive-definite measure e~ x2 (T-L m (x)) 2 dx, is an EOP 
system X m of codimension m. As shown in appendix A, these polynomials satisfy the 
second-order differential equation 



6 



3 Harmonic oscillator rational extensions in second- 
order SUSYQM 



In second-order SUSYQM (see, e.g., [27J S6j H7J SHJ HHJ [50] and references quoted therein), 
the first-order differential operators A\ A of equation (12.11) are replaced by second-order 
ones 

d 2 d d 2 d 

A j = — -2p(x)— + q(x), A = — + 2p(x) — + 2p'{x)+q{x), (3.1) 
dx z dx dx z dx 

which intertwine with two partner Hamiltonians 

= ~ + W(x), = + V&(x), (3.2) 

dx z dx A 

as AH^ = H^A and A^H^ 2 ' = H^A^. As a consequence, the functions p(x), q(x) and 

the potentials V^ l,2 \x) are constrained by the relations 

// / / \ 2 2 
2 P , ( P\ c 



q(x) = —p +p - — + I — 



2p \2pJ lQp 2 

// / / \ 2 2 ^ ' 

p p \ c 



16p 



2' 



where c is some integration constant. 

In the reducible case (corresponding to c 6 K) that we consider here, the operators A* 
and A can be factorized into products of first-order differential operators A^ = —d/dx + 
W®(x), A« = d/dx + W®(x), i = l,2, as A* = A^UW and ^ = y4 (2) A (1) , while 
and AA ] become .4U = (#(i) + |) (#« _ |), AA* = (# (2) + f ) (# (2) - f) , respectively. 
The first set of operators A^ , A 1 - 1 ) can be associated with partner Hamiltonians of the 
same type as those given in (12.11) . namely = A^A^ 1 ' and H^~' = A^ A^\ with 

corresponding potentials V^ + \x), V^'(x) and a factorization energy E\. In the same way, 
for the second set of operators A^\ A^ 2 \ we may consider = A^ A^ 2 \ H(~) = 

A^ A^\ with respective potentials V^ + '(x), V^~'{x) and a factorization energy E^. 

The relation between both approaches is given by = + § and H & = #R _ | 
with an intermediate Hamiltonian H = —d 2 /dx 2 + V(x) = + § = — f , while the 
constant c is related to the two factorization energies through c = E\ — E 2 and the function 
p(x) can be expressed in terms of the two superpotentials as p(x) = (W^ + W^)/2. 



These superpotentials can be obtained from seed eigenfunctions (f>^(x) and cf)^ 2 \x) of 



ff« and as w®\ 



1, 2, respectively. This is equivalent to 



considering two seed eigenfunctions (f>i(x), (t>%{x) of the starting Hamiltonian such 
that = 0i (x) and (2) (x) = A^^x) = W(0i, 2 )/0i, where W(0i,0 2 ) denotes the 

Wronskian of 4>i(x) and 2 In terms of the latter, one can write 



V i2 \x) = V^(x) - 2^1ogW(0 1 ,0 2 ). 



(3.4) 



To be an acceptable quantum potential, V^(x) must be nonsingular in the domain of 
definition of V^'(x), which implies that the Wronskian W(0i,02) must be nodeless. 

In the harmonic oscillator case, let us start with a pair of partner potentials V^(x) 
given by equations (12.41) and (12.81) . where we replace m by m\ (and similarly in all equations 
of section 2). The first factorization energy E\ = —1m\ — 1 being less than the ground-state 
one, a nodeless Wronskian W(0i, 2 ) can be obtained by assuming that E 2 < Ex and that 
2 (x) has a single zero on R (see, e.g., [27J 113 SB])- This can be achieved here by taking 
E 2 = — 2m 2 — 1 and 2 (^) = (pm 2 ( x ) — Hm 2 (. x ) ex2 ^ 2 with m 2 odd and such that m 2 > m\. 
The Wronskian therefore becomes 



W(0i,0 2 ) = W(0 mi) m2 ) = ^(x)e x ' 



(3.5) 



where 



~~ W("H mi ,"H m2 ) — / Hm 1 y H m2 ~ ^ mi ^m 2 — 2(m 2 ^ mi ^m 2 -i — m{H mi _i'H m2 ) (3.6) 

is a /ith-degree polynomial with /x = mi +m 2 — 1, its highest-degree term being 2 M+1 (m 2 — 
niijx^ . 

As a result, we can write the two superpotentials as 



a// 

' *"Tf. 



X) = —x — 



W^ 2 \ 



x) 



—x + 



while the potentials obtained in the two equivalent approaches are given by 
V M( x ) = x 2 , 



(3.7) 



V { -\x) = V {+ \x) = x 2 -2 



mi 



V { -\x) = x 1 



9v 



nil i 

"■mi 



7)11 

2" 

' -4. 



2- 



(3i 



9m 



8 



and 



x 



X + mi + 777,2 + 1) 



X 



V{x) 
V^(x)=x 2 



nil I 



mi 



H 



mi 



777,1 + 7772 



+ 7771 + 7772 - 3, 



(3.9) 



respectively. 

As examples of harmonic oscillator rational extensions obtained in second-order 
SUSYQM, let us quote 

32a; 2 384a; 2 



V (2) ( 



y(2) ( 



X) 



x 2 + 



x 2 + 



4a; 4 + 3 (Ax 4 + 3) 2 
24(4x 4 + 5) 



+ 2 if 777-1 = 2, 7772 = 3, 

160(28a; 4 + 20a; 2 + 5) 



8a; 6 + 20a; 4 + 10a; 2 + 5 (8x 6 + 20a; 4 + 10a; 2 + 5) 2 

if 777i = 2, 7772 = 5, 

16(16x 6 + 28a; 4 + 140a; 2 - 749) 



+ 4 



(3.10) 



(3.11) 



V (2 \ 



x 



x 2 + 



16a; 8 + 112a; 6 + 168x 4 + 84a; 2 + 21 
896(1072x 6 + 1932a; 4 + 1008a; 2 + 273) 



V <2) { 



(16a; 8 + 112a; 6 + 168a; 4 + 84x 2 + 21) 2 
64(4a; 6 + 4a; 4 - 13a; 2 + 112) 



+ 6 if mi = 2, m 2 = 7, (3.12) 



x 2 4 



16a; 8 + 64a; 6 + 120a; 4 + 45 
1024(328a; 6 + 1020a; 4 + 90a; 2 + 315) 
(16a; 8 + 64a; 6 + 120a; 4 + 45) 2 



6 if 777i = 4, 7772 = 5, 



(3.13) 



x 2 + 



8 (80a; 8 + 272x 6 + 352x 4 + 284a; 2 + 239) 



32a; 10 + 272a; 8 + 784a; 6 + 840a; 4 + 210x 2 + 105 
64(13488x 8 + 68992a; 6 + 103320x 4 + 40320a; 2 + 4515) 



(32a; 10 + 272a; 8 + 784a; 6 + 840a; 4 + 210x 2 + 105) 2 

if 777i = 4, 7772 = 7. 



+ 



(3.14) 



It is worth observing here that some of these potentials have already made their occurrence 
in another context [28l l4"2l l5Tj from rational solutions of the Painleve IV equation [52]. The 
first one of them, for instance, is a member of the shifted harmonic oscillator hyper confluent 
hierarchy of [12] and also appears in combination with a standard oscillator in the two- 
dimensional superintegrable system studied in example B of [5T] . 
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From SUSYQM, we directly get the energy spectra of H and in the form 
= 2v + mi + m 2 + 2, v = 0, 1, 2, . . . , 
S„ = 2v + mi+m 2 + 2, v = -mi - 1, 0, 1, 2, . . . , (3.15) 
£f, 2) = 2z/ + mi + m 2 + 2, z/ = -m 2 - 1, -m 1 - 1, 0, 1, 2, ... . 

In particular, we note that, as expected, the two levels of energy E_ m2 _ x = m\ — m 2 and 

(2) 

E_'_ x = m 2 -m 1 , have been added below the usual harmonic oscillator spectrum starting 

(2) 

with Eq = mi + ^2 + 2. 

The corresponding wavefunctions ip^\x) and ip v {x) being given by the right-hand sides 
of equations (I2.3P and (I2.14p (with m replaced by mi), respectively, it only remains to deter- 
mine those of . Apart from the ground-state one, which is proportional to (<j)^(x)) , 
they can be most easily found by acting with the operator on ip u (x) (see appendix A). 
The results read 

--X 2 

^\x)=^^—y^(x), n = z/+u+2, v = -m 2 -l, -mi-1, 0, 1, 2, . . . , (3.16) 
where yn(x) is an rath-degree polynomial in x, defined by 

yml+m 2 +v+i( x ) = ( m 2 - m 1 )'H mi 'Hm 2 H u+ i + 2[mi(m 2 + v + l)H mi -iH 
- m 2 (mi + v + l)l-L mi l-L m2 -i}H u , v = 0, 1, 2, . 
and the normalization coefficient is given by 



(3.17) 



)ll 2 



M {2) - ( 

- /V -m 2 -l — I 



2 m2+1 m 2 !(m 2 — mi) 



1/2 



2 mi-l mi | \l/2 (3.18) 



'7T 

\ I ' 

AA 2) = [0^ +2 (i/ + mi + l){v + m 2 + I)*/!]" 1 / 2 , i/ = 0, 1, 2, . . . . 

The nth-degree polynomials y^\x), n = mi, m 2 , mi+m 2 -|-l, mi + m 2 + 2, . . . , form an 
orthogonal and complete set with respect to the positive-definite measure e~ x2 (g^x)) 2 dx. 
We therefore conclude that such a set is a new EOP system X mi>m2 of codimension u = 

10 



rrii + m 2 — 1- The second-order differential equation satisfied by yn(x) can be written as 

(for details see appendix A) 

d 2 



dx 2 



-21 x 
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W dx J 



0. 



(3.19) 



where g^_ 2 (x) is a (/i — 2)th-degree polynomial defined by 

9»-2{x) = W{U' mV U' m2 ). (3.20) 

As a last point in this section, it is interesting to note that the order of the seed 
functions <pi and <p 2 may be changed without affecting the final results, only intermediate 
potentials and Hamiltonians being modified. So, on taking <f>i(x) = <f>2(x) = T-i m2 {x)e x2 
[m 2 odd) and 02 (^) = 4>ii. x ) — 1~Lmi( x ) ex2 ( m i even and such that rri\ < m 2 ), we obtain 
A« = d/dx + W {i \x), W^(x) = -(log0 (i) (x))', z = 1, 2, with 0W(x) = Mx) = <j) 2 {x) 
and <j>P>(x) = AWfaix) = W(<fo,<f>i)/<fo. Equations (JHU) and (J33D now become 



V {+ \x) = V {+ \x) 
V { ~\x) = V {+ \x) 



x 



x 2 - 2 



V^~\x) = V { -\x) = x 2 -2 



nil i 



0£ 



7«2 



m.2 
/ \ 2' 



(3.21) 



it 
9,, 



and 



^/ (1) (:r) = V (1) (x) = x 2 + m 1 + m 2 + 1, 



^(x) = x 1 



nil/ 



1112 



rri2 



u 



m 2 



+ mi + m 2 - 1, 



(3.22) 



V'( 2 )( a ;) = y®(x) = x 2 



+ mi + m 2 - 3, 



% _ 
respectively. 

The resulting alternative factorizations A = A^'A- 1 ' = A^A^ can be summarized in 
the following commutative diagram 



aw 



H 



AC 2 ) 



(3.23) 



Although only formally defined (since V(x) is singular at x — 0), the Hamiltonian if will 
prove very useful in the next section for constructing a new set of ladder operators for H ^ . 
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4 Ladder operators for harmonic oscillator rational 
extensions 

In the present section, we will start by reviewing the standard way of constructing ladder 
operators for in SUSYQM before introducing a new one and studying its properties. 

4.1 Standard ladder operators for 

The usual procedure for building ladder operators for a SUSYQM partner Hamiltonian 
consists in combining those of the starting Hamiltonian with the supercharge operators |26j . 
In the present case, from the standard harmonic oscillator raising and lowering operators 

a) = — — + x, a = — + x, [a, aM = 2, (4.1) 
ax ax 

valid for #«, those of are obtained in the form 

b ] = Aa ] A\ b = AaA* (4.2) 

and are therefore fifth-order differential operators. 

The operators H^ 2 \ and b satisfy the commutation relations 

[# (2) , 6 f ] = 2b\ [H {2 \ b] = -26, [b, 6 f ] = P(# (2) + 2) - P(# (2) ), (4.3) 

where P(H^) is a fifth-order polynomial in H^ 2 \ which can be factorized as 

P(H {2) ) = (H {2) -m x -m 2 - 2)(H {2) -m 1 +m 2 - 2){H {2) -m 1 + m 2 ) 

x (H {2) +m 1 -m 2 - 2)(H {2) + m x - m 2 ). (4.4) 

Hence they close a polynomial Heisenberg algebra (PHA) of fourth order. 

Some information on the unitary irreducible representations of this PHA can be obtained 
by determining the zero modes of b and b^ (i.e., the states that satisfy either hip = or 
tfijj = 0). Those of b can be deduced from the vanishing of the norm of In/j, which involves 
the average of the operator product tfb = P(H^ 2 '). Their energies are therefore given by 
m i — m 2, m i — m 2 + 2, m 2 — mi, m 2 — m 1 + 2 and mi + m 2 + 2. From (13.151) . we however 
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know that a physical state can be associated with only three of these energies, namely 



E 



(2) 

-OT2 — 1 



mi — m 2 , E_ _ x = m 2 — vti\ and E^' = mi + m 2 + 2. In the same way, the zero 



modes of 6+ can be obtained from the average of 66+ = P(# (2) + 2), leading to the energies 

(2) 

m i — m 2 — 2, mi — m 2 , m 2 — mi — 2, m 2 — mi and mi + m 2 . Here only E_'_ x = mi — m 2 

(2) 

and -Ei^-i = m 2 — mi correspond to a physical state. We conclude that the PHA, defined 
in (14. 3p and (I4.4p . has two one-dimensional unitary irreducible representations spanned by 



the singlets {ip^^^} and {^m^J, respectively, and a single infinite-dimensional one 

(2) 

spanned by {ipl \ v — 0, 1, 2, . . .}. 

4.2 New ladder operators for 

The construction of alternative ladder operators for is based upon the possibility of 
going from the first intermediate Hamiltonian H of (13.231) to the second (formal) one H (up 
to some additive constant) by a chain of I = m 2 — mi first-order SUSYQM transformations 
characterized by the supercharges 



■,(2) 



H! 

Wi(x) = x + 



mi+i— 1 



a = ^- + mix) 

ax 



mi+i 



(4.5) 



On defining 
Hi = - 



dx 2 



+ x z 



njii 

rl mi+i-l 



1,2, 



H 



mi+i— 1 



nil 

^mi +i— 1 



1,2,. ..,£+1, 



(4.6) 



it is indeed straightforward to show that 



A\k\ = H u AiA\ = H i+1 + 2, 2 = 1,2,... 



(4.7) 



which implies that Hi and Hi +X + 2 intertwine with Ai and A] as AiHi = (-Hj+i + 2)v4j and 
^ij = A](H i+ i + 2) for i = 1, 2, . . . , t Since 



H = H 1 + m 1 +m 2 + 2, H = H l+X + mi + m 2 + 2 



(4.i 



we infer that 



i^--i 2 ii#= (i? + 2£)i^--i 2 ii, Fi t 1 i^---iJ = i t 1 i+---i+(F + 2£), (4.9) 
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which proves the above assertion. This can be summarized in the following diagrammatic 
form 

H x \ H 2 + 2 H H 3 + 4 ^ . . . Hi + 2£ - 2 h H g+1 + 2£ (4.10) 



H Al - MM ) H + 2£ 



(4.11) 



This £th-order SUSYQM transformation from H to H + 21 can be combined with the 
two first-order ones going from to H or from H to to provide some raising and 
lowering operators for H@\ 



c t = A^A\A\ ■ ■ ■ A\A<M, c = A® M ■ ■ ■ A 2 A 1 A^, 



(4.12) 



which are (£ + 2)th-order differential operators. From the set of intertwining relations 
established here as well as in section 3, it is indeed easy to prove that 



cH V) = (#(2) + 2 £) c 



(4.13) 



or 




H (2) + 2t 



(4.14) 



Equation (14.131) may also be interpreted as a special type of (£ + 2)th-order SI, generalizing 
the third-order one considered elsewhere [5TL 1531 [54] , to which it reduces in the case where 
m-2 = mi + 1 and therefore £ — 1. 

The operators H^ 2 \ cS and c satisfy the commutation relations 

[HW,J] = 2£c\ [# (2) ,c] = -2£c, [c,c f ] = Q(H {2) + 2£) - Q(H {2) ) } (4.15) 



where 



Q(H (2) ) = {H {2) -3£) 



Y[( h(2) -2m! -1-2%) 



i=l 



(4.16) 



is a {£ + 2)th-order polynomial in H^ 2 \ Hence they close a PHA of (£ + l)th order. 

The zero modes of the annihilation operator c, deduced from the vanishing of the average 
of c^c = Q(H^), correspond to the energies —£, 2m i + £ + 2, 2m\ + £ + 4, . . . , 2m\ + 3£ and 
3£ or, in other words, mi — 7712, mi +1712 + 2, mi +m2 + 4, . . . , mi + 772,2 + 2^ and 3m2 — 3mi. 
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Only the first t + 1 of them, £ T ^ 2 _ 1 = m! — m 2 , -Eq Z; = mi + m 2 + 2, E{ } = mi + m<i + 4, 

(2) 

can be associated with some physical state. On the other hand, 
the zero modes of the creation operator c', obtained from the vanishing of the average of 
cc f = Q(if (2) + 2f), correspond to the energies -3£, 2m 1 -£ + 2, 2m 1 - £ + 4, . . . , 2m 1 + £ 
and £ or, equivalently, 3mi — 3m 2 , 3mi — m 2 + 2, 3m i — m2 + 4, . . . , mi + m2 and m2 — mi, 

(2) 

among which only the last one is a physical energy E_^ li _ l = m 2 — m\. We therefore 
obtain in this case one two-dimensional unitary irreducible representation spanned by the 

(2) (2) 

doublet of states {'0_^ 2 _ 1 , ^_ T ^ 1 _ 1 } and £ infinite-dimensional ones spanned by the states 
I J = 0, 1, 2, . . .} with i — 0, 1, . . . , £ — 1, respectively. 
This is confirmed by determining the action of the annihilation operator c on the physical 
wavefunctions ipl (x), v = — m 2 — 1, —mi — 1, 0, 1, 2, ... , given in equations (I3.16p - (l3.18p . 
This calculation, outlined in appendix B, leads to the results 



c^i m2 _i = c% = ap\ = ■ ■ ■ = ci\) K J 2 _ mx _ x = 0, 

2) 

-7712 — 1' 

(4.17) 



c^-TOi-i = (ma - mi) ( — ^i 2 i 2 _i, 

mi . i 



(2) _ i' 2 m ^ +2 u\{u + 2mi - m 2 + l)(z/ + m 2 + 1) \ 1/2 (2) 

v = m 2 — mi, m 2 — m x + 1, . . . . 
From (I4.17p . the action of the creation operator can be directly obtained by Hermitian 
conjugation and in particular we have c^-mi-i = 0- 

5 Conclusion 

In the present paper, we introduced a new EOP system by extending the type III Hermite 
X m EOP family [23] to a double-indexed one X mij?n2 , where mi is even and m 2 is odd 
and greater than mi- We showed that these new EOP can be easily expressed as linear 
combinations of mixed products of Hermite and pseudo-Hermite polynomials. Since their 
codimension \i = mi + m 2 — 1 is at least equal to four (corresponding to mi = 2 and 
m 2 = 3), such EOP do not appear in the recent classification of codimension two EOP 
systems obtainable from a classical system by a Darboux-Crum transformation [39]. This 
illustrates the interest of exploring higher codimensions too. 
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We also constructed the related ES rational extensions of the harmonic oscillator and 
observed that some of them were already known and obtained from rational solutions of 
the Painleve IV equation [2B1 1121 ED E2]- This establishes an interesting link between such 
solutions and EOP. 

Furthermore, for these rational extensions of the harmonic oscillator, we proposed a 
new set of ladder operators giving rise to a PHA of (777.2 — mi + l)th order. In contrast 
with what happens for standard ladder operators [2H1 [23 [2H1 HD 112] , the two states added 
below the harmonic oscillator spectrum belong to a single unitary irreducible representation 
of this PHA. Such an algebra also shows that the rationally-extended harmonic oscillator 
associated with the A mi m2 EOP family is endowed with a special type of (m 2 — mi+ 2)th- 
order SI, generalizing the third-order one considered in another context [511 EH 131] , to 
which it reduces whenever 7772 = 7771 + 1. 

Some interesting problems for future study are the construction of type III multi-indexed 
Hermite A mii m 2i ..., mfc EOP families and of the related harmonic oscillator extensions through 
the use of fcth-order SUSYQM, as well as the generalization of the present approach to 
type III Laguerre and Jacobi EOP families. Two-dimensional superintegrable systems with 
higher-order integrals of motion, based on the one-dimensional systems presented here, 
would also be worth investigating along the lines of two recent works [551 EE], thereby 
generalizing special case B considered in [5Tj . 

Appendix A. Wavefunctions of SUSYQM partners 

The purpose of this appendix is to provide some details on the calculation of wavefunctions 
(12121) . ([2713]) and (1336]) - (gJBD of first- and second-order SUSYQM partners, respectively. 

In such a computation, we make a repeated use of the following Hermite and pseudo- 
Hermite identities 

H n = 2nH n _ 1} H n+ i = 2xH n — 2nH n _ 1 , 

(A.l) 

2xU n + 2nH n -i, 

as well as of the differential equations 

H'; - 2xH' n + 2nH n = 0, H„ + 2xW n - 2nU n = 0. (A.2) 
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As a direct consequence of these relations, the Wronskian g^(x), defined in (I3.6p . fulfils the 
differential relations 



g'p + 2xg II = 2(m 2 - m 1 )'H mi 'H 



m 2 i 



9*1 + Zxg'n - = 2g 



H-2, 



(A.3) 



where g M _ 2 (x), defined in (13.201) . can also be written as g M _ 2 = —2mxH mi 'H' m2 + 
2m 2 % mi 'H m2 . 

For the first-order SUSYQM partner, the excited-state wavefunctions (I2.12p can be 
derived from the SUSYQM property 



H, 



x) 



X 



*tfr> 



-±ai 2 

e z 

% m {x) 



d 

dx 



-2x 



— 7~Li 



d _ 



#„(:r)e 2 



(A.4) 



which directly leads to the desired result after using (1A.1I) . On the other hand, the nor- 
malized ground-state wavefunction (I2.13P has been taken from [23] • It is worth noting here 
that in the same work (23], the EOP y^+ m+1 (x), v — 0, 1, 2, . . . , given in (|2.12p . have 



0, 1, 2, . 

been instead written as a linear combination of m + 1 Hermite polynomials. Finally, the 
second-order differential equation ( 12. 16ft . satisfied by yn(x), can be easily obtained by 
inserting equation (12.141) in the Schrodinger equation for V^~\x) and using (1A.2I) . 

For the second-order SUSYQM partner, the ground-state wavefunction ip_m 2 -i( x ) comes 
from the inverse of (j)^(x), leading to the expression of yml(x) in (13.171) . On the other hand, 
the excited-state wavefunctions result from the action of (with W^ 2 \x) given in (13. 7p ) 
on the intermediate Hamiltonian wavefunctions ip v (x), v = —mi — 1, 0, 1, 2, ... . On taking 
the first relation in ( 1A.3I) into account, this yields 



A {2) i) u {x) oc 



e 2 X 1 



g^-2(m 2 - m 1 )'H mi 'H m2 



(mi) 
»i/+mi+l 



(A.5) 



from which the expression of yml{x) (corresponding to v — —m\ — 1) in (13.171) follows. 
For higher values of v [y = 0, 1, 2, . . . ), we have to insert definition (13. 6 p of g^{x) in the 
right-hand side of equation (1A.5I) . Among the resulting three terms inside square brackets, 
the first and third ones contain "H mi as a factor, while the second one does not. After using 



17 



the second relation in (I2.12p as well as equation (lA.lj) . this second term can however be 
rewritten as 

- 2miH mi -iH m2 ^y^lll 1+1 = Am x {m x + v + l) / H mi 'H mx -iH m2 H l ,. (A.6) 

Hence all three terms in the numerator of ( 1A.5I) contain the factor 1-L mi , which cancels 
the same in the denominator. On employing (1A.1I) again, the remaining expression can be 
further transformed into the third relation in equation f)3.17p . 

The normalization coefficient v = —mi — 1, 0, 1, 2, . . . , of the excited-state 

wavefunctions can be deduced from that of the intermediate Hamiltonian wavefunctions, 
from which such wavefunctions are derived, by multiplying the latter by the extra factor 
(-E^) 1//2 = [2(z/ + m 2 + 1)]~ 1//2 . This yields the second and third relations in f)3.18p . while 

(2) (2) 

the first one in the same is obtained from the second by observing that ip_ m2 _± and ip^'i 
have a similar form up to a permutation of mi and m 2 . 

Finally, the differential equation ( 13 .191) comes from inserting equation ( 13. 161) in the 
Schrodinger equation for V^(x) and applying the second relation in (1A.3[) . 

Appendix B. Action of the new annihilation operator c 

The purpose of this appendix is to prove that the action of the operator c, defined in (14. 12ft . 
on the wavefunctions iph (x), v = —m 2 — 1, —m 1 — 1, 0, 1, 2, ... , given in (I3.16p - fl3.18p . is 
provided by equation f )4.17p . Such a result actually agrees with SUSYQM predictions, but 
since the operators Ai, i — 1, 2, . . . , £, and A^ 2 \ used in definition (I4.12p . are only formally 
defined, we find it appropriate to make such an explicit check. 

Since is a well-defined operator on M, we may use its SUSYQM properties yielding 

A^ = {\ if, = -m 2 -l, 

]y 2 + ™2 + lM, if v = -m x - 1, 0, 1, 2, .... 

This shows that c^i 2 ^ 2 „ 1 = and 

--X 2 

cx A^Ae ■ ■ ■ A 2 A l \^y i Jlll i+1 , v = - mi - 1, 0, 1, 2, . . . . (B.2) 



18 



To calculate the action of the product Ag - • • A 2 Ai in flB.2j) . we repeatedly employ the 
relation 



At 



e 2 



14 



(mi+i— 1) 
'Vu+mi+l 



-±x 2 

e * 



u mi 



+i 



Hm 1 +i— 7 2(mi + i)l-t mi+ i-i 

ax 



(mi+i— 1) 
Vv+mi+l 



- z + 1)— ^ +mi j 1; i = l,2 



which results from (12.121) . (14. 5p . (lA.lj) and the use of the intermediate equation 
d 



(mi+i-l) 

dx Jv+mi+1 



-.'/;,+;,, , +i" ' = -2(v + mi+ 1 !'#,,,. ., .///, , . i* i = 1, 2, 
We finally arrive at 



e# oc A® 



-1/ fv+mi+l j 14 a 14 ^f+mi+lJ 



z/ = - mi - 1,0, 1,2, ... , 



which can be easily shown to be given by 



-It 2 

,(2) „ e 5 „» 



; OC V = ~ m l ~ X ' °' X > 2 ' • • • 

57* 



on using f l2TT2|) . flATT]) and ([Q]) . 

After collecting all coefficients, the result reads 

' 2£[(2m x + 2) (2m! + 4) ■ • • (2m! + 2£)] 1 / 2 ^ ) m2 _ 1 if u = -m x - 1, 
if z/ = 0,l,...,£-l, 

[2i/(2i/ - 2) • • • {2v - 21 + 2){2u + 2rn 1 - 2t + 2)] 1 / 2 

x (2z/ + 2mi + 2£ + 2) 1 /2^ if z/ = £, t + 1, £ + 2, . . ., 

which is equivalent to equation (14. 17ft when taking the definition of i into account. 



c^ 2) = < 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



19 



References 

[1] Gomez-Ullate D, Kamran N and Milson R 2009 J. Math. Anal. Appl. 359 352 

[2] Gomez-Ullate D, Kamran N and Milson R 2010 J. Approx. Theory 162 987 

[3] Quesne C 2008 J. Phys. A: Math. Theor. 41 392001 

[4] Gendenshtein L E 1983 JETP Lett. 38 356 

[5] Cooper F, Khare A and Sukhatme U 2000 Supersymmetry in Quantum Mechanics 
(Singapore: World Scientific) 

[6] Gomez-Ullate D, Kamran N and Milson R 2004 J. Phys. A: Math. Gen. 37 1789 

[7] Gomez-Ullate D, Kamran N and Milson R 2004 J. Phys. A: Math. Gen. 37 10065 

[8] Bagchi B, Quesne C and Roychoudhury R 2009 Pramana J. Phys. 73 337 

[9] Quesne C 2009 SIGMA 5 084 

[10] Odake S and Sasaki R 2009 Phys. Lett. B 679 414 

[11] Odake S and Sasaki R 2010 Phys. Lett. B 684 173 

[12] Odake S and Sasaki R 2010 J. Math. Phys. 51 053513 

[13] Ho C-L, Odake S and Sasaki R 2011 SIGMA 7 107 

[14] Sasaki R, Tsujimoto S and Zhedanov A 2010 J. Phys. A: Math. Theor. 43 315204 

[15] Gomez-Ullate D, Kamran N and Milson R 2010 J. Phys. A: Math. Theor. 43 434016 

[16] Gomez-Ullate D, Kamran N and Milson R 2012 Contemp. Math. 563 51 

[17] Grandati Y 2011 Ann. Phys., NY 326 2074 

[18] Ho C-L 2011 Prog. Theor. Phys. 126 185 

[19] Dubov S Y, Eleonskii V M and Kulagin N E 1992 Sov. Phys. JETP 75 446 

20 



[20] Dubov S Y, Eleonskii V M and Kulagin N E 1994 Chaos 4 47 

[21] Junker G and Roy P 1997 Phys. Lett. A 232 155 

[22] Carinena J F, Perelomov A M, Ranada M F and Santander M 2008 J. Phys. A: Math. 
Theor. 41 085301 

[23] Fellows J M and Smith R A 2009 J. Phys. A: Math. Theor. 42 335303 

[24] Sukumar C V 1985 J. Phys. A: Math. Gen. 18 2917 

[25] Junker G and Roy P 1998 Ann. Phys., NY 270 155 

[26] Fernandez C D J and Hussin V 1999 J. Phys. A: Math. Gen. 32 3603 

[27] Fernandez C D J and Fernandez-Garcia N 2005 AIP Conf. Proc. vol 744 (Melville, 
NY: Amer. Inst. Phys.) p 236 

[28] Bermiidez D and Fernandez C D J 2011 SIGMA 7 025 

[29] Grandati Y 2011 J. Math. Phys. 52 103505 

[30] Ho C-L 2011 J. Math. Phys. 52 122107 

[31] Gomez-Ullate D, Kamran N and Milson R 2012 J. Math. Anal. Appl. 387 410 

[32] Odake S and Sasaki R 2011 Phys. Lett. B 702 164 

[33] Quesne C 2011 Mod. Phys. Lett. A 26 1843 

[34] Quesne C 2011 Int. J. Mod. Phys. A 26 5337 

[35] Grandati Y 2012 Ann. Phys., NY 327 2411 

[36] Quesne C 2012 Int. J. Mod. Phys. A 27 1250073 

[37] Quesne C 2012 SIGMA 8 080 

[38] Grandati Y 2012 Phys. Lett. A 376 2866 

21 



[39] Gomez-Ullate D, Kamran N and Milson R 2012 A conjecture on exceptional orthogonal 
polynomials larXiv: 1 203.6857 

[40] Sasaki R and Takemura K 2012 SIGMA 8 085 

[41] Carballo J M, Fernandez C D J, Negro J and Nieto L M 2004 J. Phys. A: Math. Gen. 
37 10349 

[42] Mateo J and Negro J 2008 J. Phys. A: Math. Theor. 41 045204 

[43] Abramowitz M and Stegun I A 1965 Handbook of Mathematical Functions (New York: 
Dover) 

[44] Noumi M and Yamada Y 1999 Nagoya Math. J. 153 53 
[45] Clarkson P A 2003 J. Math. Phys. 44 5350 

[46] Andrianov A A, Ioffe M V and Nishnianidze D N 1995 Phys. Lett. A 201 103 
[47] Bagrov V G and Samsonov B F 1995 Theor. Math. Phys. 104 1051 
[48] Samsonov B F 1996 Mod. Phys. Lett. A 11 1563 

[49] Bagchi B, Ganguly A, Bhaumik D and Mitra A 1999 Mod. Phys. Lett. A 14 27 
[50] Aoyama H, Sato M and Tanaka T 2001 Nucl. Phys. B 619 105 
[51] Marquette I 2009 J. Math. Phys. 50 095202 

[52] Gromak V, Laine I and Shimomura S 2002 Painleve Differential Equations in the 
Complex Plane (Berlin: Walter de Gruyter) 

[53] Andrianov A, Cannata F, Ioffe M and Nishnianidze D 2000 Phys. Lett. A 266 341 

[54] Ioffe M V and Nishnianidze D N 2004 Phys. Lett. A 327 425 

[55] Post S, Tsujimoto S and Vinet L 2012 J. Phys. A: Math. Theor. 45 405202 

[56] Marquette I and Quesne C 2012 New families of superintegrable systems from Hermite 
and Laguerre exceptional orthogonal polynomials ar Xiv:1211 .2957 

22 



